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We perform a global fit to the CKM unitarity triangle using the latest experimental and theo- 
retical constraints. We present results for three different sets of constraints: the standard inputs 
used by CKMfitter [1] and UTFit [2], the standard inputs minus \Vub\, and the standard inputs 
minus both \V„b\ and \Vcb\ [3]. For the required nonperturbative weak matrix elements, we use 
three-flavor lattice QCD averages from www . latticeaverages . org; these have been up- 
dated from Ref. [4] to reflect all available lattice calculations as of the ^'End of 2010". Given 
current theoretical and experimental inputs, we observe an approximately 3(J tension in the CKM 
unitarity triangle that can be interpreted as sign of physics beyond the standard model in the 
flavor sector. Using a model-independent parameterization of new physics effects, we test the 
compatibility of new physics in kaon mixing, in B-mixing, or in B ^ TV decay with the current 
data. Although the tension could be accommodated with each hypothesis, the scenarios with new 
physics in B-mixing or, to a lesser extent, in B ^ TV decay are strongly preferred. Finally, we 
interpret these results in terms contributions to AS' = 2 and AS ~ 2 four-fermion operators. We 
find that the preferred scale of new physics (with Standard Model like couplings) is in the few 
hundred GeV range. 
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1. Motivation 

The B-factories and the Tevatron have produced a remarkable wealth of data needed to deter- 
mine elements of the Cabibbo-Kobayashi-Maskawa (CKM) matrix and to search for new physics 
beyond the Standard Model CKM framework. Despite the great experimental success of the Stan- 
dard Model, there is now considerable evidence for physics beyond the Standard Model, such as 
dark matter, dark energy, and neutrino masses. Generic new physics models to explain such phe- 
nomena also lead to additional CP- violating phases beyond the single one in the Standard Model; 
this would lead to apparent inconsistencies between independent determinations of the CKM ma- 
trix elements. Although there is presently reasonably good agreement with the Standard Model 
prediction of a single CP-violating phase, as measured by global fits of the CKM unitarity triangle, 
some tensions have been observed [5-11]. In this work we use the latest theoretical and experi- 
mental inputs to quantify the tension with the Standard Model via a global fit to the CKM unitarity 
triangle and then identify within a largely model-independent framework the most likely sources 
of the new physics. 

2. Unitarity Triangle Fit Preliminaries 
2.1 Inputs 

The standard analysis of the unitarity triangle involves a simultaneous fit to several quantities: 
Ek, ^^Bd' AMg,, time-dependent CP asymmetry in B — )• J/^Ks iSif,K = sin(2jS), where j8 is the 
phase of V*^),^ direct CP asymmetries in B — )■ d'*'^'*) (7 is the phase of V,*^) time dependent 
CP asymmetries in B — ;> {nn,pp,pn) (a = tt - jS - 7), BR(B tv), \Vub\ and \Vch\ (from both 
inclusive and exclusive — )• {u,c)lv with 1 = e,v). We summarize the relevant inputs required for 
this analysis in Table 1 . 

For the nonperturbative weak matrix elements we use averages of three-flavor lattice QCD cal- 
culations from Ref . [4] updated to reflect all results documented in proceedings or publications as 
of the ''End of 2010" . We do not include two-flavor lattice calculations in our averages because of 
the unknown systematic en^or due to neglecting dynamical strange quark effects. We treat both sta- 
tistical and systematic en^ors as following a Gaussian distribution. In performing the averages, we 
account for correlations between different lattice calculations of the same quantity in a reasonable 
but conservative manner; whenever there is reason to believe that an error is correlated between 
two results, we assume that the degree of correlation is 100%. We obtain errors in pion and kaon 
matrix elements that are consistent with those of the Flavianet Lattice Averaging Group when we 
use the same inputs [20], despite the fact that they use a different method for combining systematic 
errors between lattice calculations. FLAG, however, has not yet updated their averages to reflect 
recent results from Lattice 2010 nor do they present averages of B- or D- meson quantities. 

There are several choices for how to implement the constraints from B ^ TV leptonic decay 
and Bd s mixing (AM5, and AM^,.,) on the unitarity triangle because one can parameterize the 
nonperturbative weak matrix element contributions to these quantities in different ways. Certain 
combinations of lattice inputs are preferable, however, because they minimize correlations between 

'a discussion of penguin pollution in S^/k can be found in Ref. [12]; see also Refs. [13-19]. 
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the three different unitarity triangle constraints so that they can safely be neglected in the global 
fit. Let us now summarize the main considerations that lead to a reasonable choice of uncorrelated 
inputs. First of all it is important to include only one input with mass dimension 1 in order to 
eliminate correlations due to the determination of the lattice scale. Another important consideration 
is that the largest source of uncertainty in the SU (3)-breaking ratios, and fsj Ib^, is the chiral 
extrapolation. Because the chiral logarithms are larger when the quark masses are lighter, the chiral 
extrapolation in the SU (3) -breaking ratios is primarily due to the chiral extrapolation in the 
quantities (i.e. /g, and B^). These ratios are, therefore, more correlated with rather than with B^ 
quantities. Finally we note that the decay constants /b, and fs^ have a stronger chiral extrapolation 
than the B-parameters B^ and B^. In view of these considerations we choose to describe Bq mixing 
in terms of fB,\fB^s and t,; once this choice is made, the additional input required to describe 
B — > TV has to be B^ (because /g^ has mass dimension and is somewhat correlated with t,). For 
completeness we point out that there is an alternative choice of inputs {fsj fB^, B^/B^, /b, and B^) 
for which correlations are again fairly small. 



In our analysis, we write the unitarity triangle constraints in terms of the lattice inputs fs, v B 
t, , and B^ such that the the unitarity triangle constraints are 



An additional advantage to this choice of inputs is that it allows us to use all existing experimental 
data to obtain a prediction for /bj from all other constraints on the unitarity triangle that can be 
compared to the direct lattice calculation. 

The determinations of \Vcb\ and \Vub\ from inclusive and exclusive decays are problematic 
because they differ at the 2.1 and 1.7 a levels, respectively (note that if we remove the additional 
10% model uncertainty on inclusive \Vuh\, the discrepancy rises to the 3.3 a level). 

The determinations of \Vch\ and \Vub\ from inclusive and exclusive decays are problematic 
because they differ at the 2.1 and 1.7 s levels, respectively (note that if we remove the additional 
10% model uncertainty on inclusive \Vuh\, the discrepancy rises to the 3.3a level). For this reason, 
when we combine the inclusive and exclusive determinations of \Vcb\ and \Vuh\ we inflate the errors 
by the square root of the chi-square per degree-of-freedom (as prescribed by the PDG); the resulting 
averages are given in the bottom panel of Table 1. Furthermore, in addition to the standard fit in 
which all measurements are included, we consider two additional scenarios in which we remove 
\Vuh\ and \Vch\ from the chi-square. The strategy for removing \Vch\ by combining the constraints 
from Sk, AMb,, and BR{B — > Tv) was proposed and is described in detail in Ref. [3]. 

2.2 Interpretation as New Physics 

We interpret the observed tensions in the global unitarity triangle fit as contributions from new 
physics in either in kaon mixing, S^-mixing, or B — )• tv; details of the analysis method are given 
in Ref. [11]. We adopt a model-independent parametrization of new physics effects in the three 





(2.1) 
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lK:folexcl = (39.5 ± 1.0) X 10-3 

Bjf = 0.737 ±0.020 
/b = (205 ± 12) MeV 
Bb, = 1.26±0.11 

^=fB,VB,/{fB,VBd) = 1.237 ±0.032 



|Ki,lexci = (3-12±0.26)xlO-3 

K-e = 0.94 ±0.02 

/b, = (250±12)MeV 

% = 1.33 ±0.06 

/B,V^=(288±15)MeV 

/By/B, = 1.215 ±0.019 



|V,^|.^^j = (41 .68 ± 0.44 ± 0.09 ± 0.58) x 10"^ [21] a = (89.5 ± 4.3)° 
\VubLa = (4.34 ±0.16^^-15 ±0.43) x lO'^ [21] 771 = 1.51 ±0.24 [22] 
BR(S^ TV) = (1.68±0.31) X 10-^ [23-25] 



AniB, = (0.507 ±0.005) ps'^ [27] 
Ams^ = ( 17.77 ±0. 10 ± 0.07) ps-i [29] 
m,,p„,, = (172.4±1.2)GeV[31] 
mdmc) = (1.268 ±0.009) GeV [34] 
Ek = (2.229 ±0.012) x lO^^ [36] 



5y,jf, = 0.668 ±0.023 [26] 
y=(78±12)° [2,28] 

7]2 = 0.5765 ±0.0065 [30] 

T73 =0.494 ±0.046 [32,33] 

TjB = 0.551 ±0.007 [35] 

A = 0.2253 ±0.0009 [37] 



l^ci,|avg = (40.77±0.81)xl0-3 



|Kfelavg = (3.37±0.49)xl0-3 



Table 1: Lattice QCD and other inputs to the unitarity triangle analysis. The determination of a is obtained 
from a combined isospin analysis of B ^ i'^'^, pp, p7f) branching ratios and CP asymmetries [27]. Refer- 
ences for the lattice-QCD results entering the averages in the the upper panel can be found in Ref. [4] with 
updates at www . latticeaverages . org. Updated lattice averages for other quantities that do not enter 
the global unitarity triangle fit such as pion, kaon, and D-meson decay constants and Ught-quark masses can 
also be found at www . latticeaverages . org. 



observables: 

|eNP| = Ce lef^l , (2.2) 

A4NP^^2^2,e,^^^SM^ (2.3) 

^ ' V m^,(l±£otanj3)y ^ ' 

= rn BR{B Tv)™ , (2.5) 

where in the Standard Model (Cg, rn, r^) = 1 and = 0. In presence of non-vanishing contribu- 
tions to Srf mixing the following other observables are also affected: 

5v,jf, = sin2(/3±0rf) , (2.6) 

sin(2aeff) = sin2(a - 0j) , (2.7) 

^"' = amb:=^^' ■ ^^-^^ 
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When considering new physics in Bd mixing we allow simultaneous variations of both 6^ and r^. 
We find that new physics in \Mf2\ has a hmited effect on the tension between the direct and indirect 
determinations of sin(2j3); as a consequence, our results for and 6^ point to larger effects on the 
latter. 

Finally we interpret the constraints on the parameters Cg , r^, and in terms of generic new 
physics contributions to A5 = 2 and AB = 2 four-fermion operators. The most general effective 
Hamiltonian for B^-mixing can be written as ^ 

=^eff = ^ {VM' (t CiOi + 1 QO^i (2.9) 

where 

01 = (dLy^bL){dLy^bL) 

02 = (dRbL)(dRbL) 

o,={dibi){ibi) 

O4 = {dRbL){dLbR) 

Within the Standard Model, only the operator Oi receives a non- vanishing contribution at a high 
scale jU// ~ m,. For our analysis we assume that all new physics effects can be effectively taken 
into account by a suitable contribution to Ci : 

^.»=#(W)^cfM(^-g)o,, (2.n, 

where the minus sign has been introduced a posteriori (as we will see the the fit will point to new 
physics phases of order (p ~ 0(1)). In this parametrization A is the scale of some new physics 
model whose interactions are identical to the Standard Model with the exception of an additional 
arbitrary CP violating phase: 

Ci=cr(^l-e''P^^ . (2.12) 

When discussing new physics in the kaon sector, we will consider also a similar new physics 
contribution to the operator O4: because of RG effects and of the chiral enhancement of the matrix 
element of O4, the latter usually point to a new physics scale that is larger by a factor ^ 65 than for 
Oi case. 

3. Unitarity Triangle Fit Results and Constraints on New Physics 

In this section we present the results we obtain for the full fit and for the fits in which semilep- 
tonic decays (for the extraction of \Vub\ and Vcb\) are not used. For each set of constraints we present 
the fitted values of the CKM parameters p, ff and A. We also show the predictions for several in- 
teresting quantities (most importantly 5^,^- and BR(B TV)) that we obtain after removing the 
corresponding direct determination from the fit. Finally we interpret the observed discrepancies in 
terms of new physics in £k, S^-mixing or B ^ TV. 

■^The Hamiltonians for Sj- and A'-mixing are obtained by replacing {d, b) (i, h) and (rf, h) (rf, s), respectively. 



0\ = idRY^ibR){dR'Y^bR) 
O2 = (dLbR){dLbR) 

d, = {dlbl){d{b-^) 

O, = id",b^,){d^b-) . 



5 



Lessons for new physics from CKM studies 



Enrico Lunghi 




Figure 1: Unitarity triangle fit with all constraints included. Quantities that are not used to generate the 
black contour are grayed out. 



In Figs. 1, 3 and 5, we show the global CKM unitarity triangle fit for the three set of inputs that 
we consider (complete fit, no V„/, fit, no Vqb fit). In each figure, the black contours and p-values 
in the top, middle and bottom panels correspond to the complete fit, the fit with a new phase in 
B mixing (i.e. without using S^^/k and a) and the fit with new physics in B — ^ TV (i.e. without 
using BR(B — > Tv)), respectively.^ In the fits with a new phase in B mixing we also show the fit 
predictions for sin(2j8) and /b; in the fits with new physics in B — > TV, we show the fit predictions 
for BR(B — tv) and fs- Note that the individual contours in Figs. 1-5 never use the same input 

■'Note that in the no Vqi, fit we define the new physics in S — !• TV scenario by removing B — !• TV, AMb, and Ek from 
the fit. 
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Figure 2: Unitarity triangle fit with all constraints included: new physics analysis. 



twice in order to minimize the correlations between constraints: in particular, the 5 — >• TV allowed 
area is obtained by using AM^^ instead of the direct determination of \Vcb\- 

In Figs. 2, 4 and 6 we show the interpretation of the tensions highlighted in the various fits 
in terms of possible new physics. In the first panel of each figure we show the result of the two- 
dimensional fit in the {6d,rd) plane and in the second panel we map this allowed region onto the 
(A,<p) plane (see Eq. (2.12)) under the assumption of new physics in 0\ only. In the third and 
fourth panels we show the allowed (A, (p) regions for the scenario with new physics in K mixing 
(the left and right panels show contributions to 0\ and O4, respectively). 

3.1 Standard Fit 

We include constraints from Ek, AM^^, AMg,, a, Sxf,K, 7, BR(S — > Tv), \Vcb\ and |V„fo|. The overall 
p- value of the Standard Model fit is p = 2.6% and the results of the fit are 

p =0.135±0.018 T7 = 0.354 ±0.013 A = 0.816±0.011 . (3.1) 

The predictions from all other information when the direct determination of the quantity is removed 
from fit are 
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\Vub\ = (3.64±0.13) X 10"^ (0.53 a) (3.2) 

5^jf = 0.795 ±0.041 (2.5 a) (3.3) 

\Vcb\ = (42.1 ±0.82) X 10"^ (1.1 a) (3.4) 

= 0.889 ±0.083 (1.9 a) (3.5) 

/i?,V^= (210.0 ± 4.3) MeV (3.6) 

BR(S ^ TV) = (0.773 ±0.096) X 10"^ (2.7 a) (3.7) 
7b, = (193.±10.)MeV (0.8 a) complete fit 



< /^^ = (198.8 ±9.9) MeV (0.4 a) without using 5^,/^ (3.8) 
/^^ = (185.9 ±8.9) MeV (1. a) without using BR(B ^ tv) 
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where we indicate the deviation from the corresponding direct determination in parentheses. The 
interpretation of the above discrepancies in terms of new physics in A'-mixing, B(/-mixing and 
B TV yields 

Ce = 1.21 ±0.12 (1.9 a, 7? = 0.056) (3.9) 
f0^ = -(4.5±2.1)'' , 

< ^ ' 2.2 a, /? = 0.13) (3.10) 

[rrf = 0.96 ±0.039 

r// = 2.22 ±0.49 (2.8 a, p = 0.25) . (3.11) 



Figure 1 summarizes these results. 



3.2 Fit without \Vub\ 

We include constraints from Ek, AM^^, AMg^, a, S^k, Y, BR(S — ^ Tv) and \Vcb\- The overall 
p-value of the Standard Model fit is p = 1 .4% and the results of the fit are 

p =0.136±0.018 T]= 0.355 ±0.013 A = 0.816±0.011 . (3.12) 
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Figure 5: Unitarity triangle fit without \Vub\ and \Vcb\- Quantities that are not used to generate the black 
contour are grayed out. 



The predictions from all other information when the direct determination of the quantity is removed 
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from fit are 



\Vub\ = (3.64±0.13) X 10"^ (0.53 a) 


(3.13) 


5v,/f = 0.861 ±0.048 (3.3 a) 


(3.14) 


iKfol = (42.04 ±0.82) X 10"^ (1.1a) 


(3.15) 


fi^f = 0.887 ±0.083 (1.9 a) 


(3.16) 


/b,-\/^= (210.2 ± 4.3) MeV 


(3.17) 


BR(B TV) = (0.778 ± 0.098) x 10""* (2.7 a) 


(3.18) 


/b, = (194.±10.)MeV (0.8 a) no fit 




fBj = (200.2 ±9.3) MeV (0.4 a) without using S^,k 


(3.19) 



/b, = (186.0 ±9.0) MeV (1.3 a) without using BR(B ^ tv) 

where we indicate the deviation from the corresponding direct determination in parentheses. The 
interpretation of the above discrepancies in terms of new physics in /^-mixing, Brf-mixing and 
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B ^ TV yields 

Ce = 1.20 ±0.12 (1.9 C7, p = 0.030) (3.20) 

<^ ^ ' 2.9 a, p = 0.30) (3.21) 

[rj = 0.97 ±0.045 

ru = l.lQ±^A'^ (2.8 a, 7? = 0.16) . (3.22) 
Figure 3 summarizes these results. 

3.3 Fit without \Vuh\ and \Vch\ 

We include constraints from e^, M^Bj, AMg^, a, S^k, 7 and BR(B — tv). The overall p-value of 
the Standard Model fit is = 1 . 1% and the results of the fit are 

p =0.139±0.018 T]= 0.349 ±0.015 A = 0.828 ±0.016 . (3.23) 

The predictions from all other information when the direct determination of the quantity is removed 
from fit are 

iKil = (3.64±0.13) X 10"^ (0.53 a) (3.24) 
S^K = 0.904 ± 0.047 (3.1 a) (3.25) 
\Vcb\ = (42.1 ±0.82) X 10"^ (1.1a) (3.26) 
= 1.11 ±0.21 (1.9 a) (3.27) 

/g^Y^= (208.2 ±4.6) MeV (3.28) 

BR(B ^ TV) = (0.763 ±0.098) X lO"'^ (2.8 a) (3.29) 
/b, = (192.±10.)MeV (1.3 a) no fit 

fs, = (201.8 ±9.3) MeV (0.2 a) without using S^^k (3-30) 
/^^ = (184.1 ±9.0) MeV (1.5 a) without using BR(B ^ tv) 

where we indicate the deviation from the corresponding direct determination in parentheses. The 
interpretation of the above discrepancies in terms of new physics in /T-mixing, B^-mixing and 
5 — )• TV yields 

Ce = 1.51 ±0.29 (2.0 a, /? = 0.027) (3.31) 
fe,, = -(10.8±3.2)° ^ 

< 3.3 a, p = 0.80 (3.32) 

[rrf = 0.95 ±0.045 

r// = 2.25 ±0.50 (2.8 a, p = 0.17) . (3.33) 
Figure 5 summarizes these results. 
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4. Discussion 

The main lessons to be discerned from the fits described in the previous sections are: 

• If we take the current experimental and theoretical errors in the inputs to the global CKM uni- 
tarity triangle fit at face value, then the CKM description of flavor and CP violation displays 
a tension at the 3 a level. 

• The use ofb^ uiv {i = e,iJ.) semileptonic decays is problematic. The (2 — 3)a discrepancy 
between the inclusive and exclusive extractions of |V„/,| when coupled to the complexity of 
theoretical methods used is cause of serious concern. Furthermore, the direct determination 
of BR(B TV) from /b and |Kfo|exci (we obtain (0.64 ±0.12) x 10""^) is 3.1a below its 
direct experimental determination (see Table 1). These two 3a tensions in the B — )• TV, /b, 
\Vub\ system could be quite naturally resolved by a shift in the direct extraction of \Vub\ due to 
either new physics (see for instance Refs. [38,39]) or further theoretical progress on inclusive 
and exclusive semileptonic decays. In the light of these considerations we believe that "fit 
without \Vub\" presented in Sec. 3.2 is the most apt to represent our present understanding of 
CP violation in the CKM. 

• The fit seems to prefer an interpretation of this 3 a tension in terms of new physics in 
mixing or B — TV rather than in Sk (see the p-values in Eqs. (3.20-3.22)). Note that even 
removing BR(S Tv) from the fit without |V„fo|, we obtain S^^ = 0.803 ±0.069 (1.9a). 
Thus the presence of new physics in B — )• TV can reduce the deviation of S^j^ from S^^^ from 
3.2a to 1.9a, but cannot completely eliminate the tension in the global unitarity triangle fit. 
In contrast, new physics in B^i mixing can completely tackle the problem: the values of p 
and T] required to eliminate the 1.9a tension, also yield a BR(B — )• tv) that is in perfect 
agreement with experiments. Therefore we conclude that the most satisfactory resolution of 
the overall tension is through new physics in Bj mixing. This conclusion is reinforced by 
the inspection of Eqs. (3.8), (3.19) and (3.30): the fit result for /g, is much closer to direct 
lattice determination in the fit with new physics in mixing (i.e. without using 5,^^) rather 
than in B — ^ TV (i.e. without using BR(S — ;> Tv))."^ 

• In terms of a new physics model whose interactions mimic closely the SM (see Eq. (2.11)), 
this tension points to a few hundred GeV mass scale. Even allowing for a generous model 
dependence in the couplings, it seems that such new particles, if the tension in the fit stands 
confirmed, cannot escape detection in direct production experiments. 
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